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I.  "MODELING  OF  IONIZED  GAS  SYSTEMS. 


p  I  C 


Program  Leader  -  Dr.  A.  V.  Phelps 


Our  studies  of  the  modeling  of  ionized  gas  systems  include  a)  the 
theoretical  prediction  of  the  conditions  fc^*  .he  growth  of  a  constriction 
in  a  weakly  ionized  gas,  such  as  found  in  a  molecular  gas  laser;  b)  the 
development  of  theories  of  radiation  transport  and  hydrodynamics  applic¬ 
able  to  a  wider  variety  of  real  situations,  such  as  expanding  and  turbu¬ 
lent  plasmas  in  geometries  other  than  the  usual  planar  geometry;  c)  the 
testing  of  radiative  transport  theories  under  controlled  laboratory  con¬ 
ditions;  and  d)  the  experimental  determination  of  the  laws  of  radiative 
scattering  and  absorption  appropriate  to  radiative  transport  problems. 

The  overall  ob.j active  of  this  program  is  to  supply  experimentally  tested 
understanding  and  computational  techniques  required  for  the  analysis  and 
prediction  of  the  role  of  a)  the  generation  and  transport  of  radiation 
in  plasmas  and  b)  the  transport  and  reactions  of  charged  particles  in  weak 
to  fully  ionized  gases .  Because  of  its  concern  with  a  number  of  simul¬ 
taneous  atomic  processes  and  with  relatively  large  numbers  of  interacting 
atoms  or  molecules ,  this  program  also  provides  an  important  link  between 
the  experimental  and  theoretical  determinations  of  cross  sections  and  rate 
coefficients  which  make  up  much  of  the  remainder  of  the  program  and  the 
technology  of  ionized  and  excited  gases  which  this  work  is  designed  to 


support. 


Stability  of  Discharges  In  Weakly  Ionized  Gases  (Dr.  E.  F.  Jaeger,  Dr.  L. 
Oster  and  Dr.  A.  V.  Phelps) 

The  objective  of  this  theoretical  project  Is  the  quantitative  Inves¬ 
tigation  of  factors  which  control  the  growth  of  Instabilities  In  weakly 
Ionized  gas  discharges,  such  as  used  in  high  power  molecular  lasers.  The 
present  phase  of  this  project  is  concerned  with  instabilities  which  result 
from  a  decrease  in  gas  density  and  the  accompanying  increase  in  the  ratio 
of  the  electric  field  to  the  gas  density  E/N  and  in  the  rate  of  ioniza¬ 
tion. 

Since  the  arrival  of  Dr.  Jaeger  in  June,  we  have  been  investigating 
the  predictions  and  validity  of  previously  published  solutions  for  the 
steady  state  spatial  distributions  of  electron  density  and  gas  temperature 
and  for  the  stability  of  these  solutions  against  fluctuations  in  discharge 
conditions.  Ue  find  agreement  with  previous  results  only  if  one  takes 
into  account  the  high  degree  of  sensitivity  of  the  results  to  changes  in 
the  material  properties,  such  as  the  ambipolar  diffusion  and  ionization 
coefficients.  An  example  of  this  sensitivity  is  shown  in  Fig.  1  where  the 
spatial  distribution  of  the  electron  density  for  a  helium  discharge  is 
plotted  as  a  function  of  radius  for  various  assumed  variations  of  the  am¬ 
bipolar  diffusion  coefficient  with  or  equivalent  variable  ET/p.  The 
value  of  n  =  1/3  is  most  nearly  correct  for  helium  in  the  usual  range  of 
E/N,  although  n  varies  significantly  among  gases. 

In  addition,  we  have  obtained  solutions  for  the  time  dependence  of 
the  electron  density,  gas  temperature,  discharge  current  and  electric 
field  for  a  range  of  applied  voltages,  source  impedances,  etc.  Some  of 
these  results  are  shown  in  Fig.  2  for  a  discharge  in  helium  with  values  ap¬ 
propriate  to  a  low  pressure,  moderate  current  discharge.  Unfortunately,  an 


examination  of  the  neutral  atom  density  profiles  for  this  case  shows  that 
the  gas  Is  expanding  with  a  velocity  which  far  exceeds  the  velocity  of 
sound  appropriate  to  the  peak  temperature.  This  calculation  confirms  our 
justification  of  this  project,  i.e.,  that  the  constant  pressure 
approximation  normally  used  lu  modeling  glow  discharges  should  not  be  used 
at  times,  pressures  and  dimensions  of  interest  in  many  laser  problems.  We 
are  currently  formulating  the  hydrodynamic  and  charged  particle  equatlans 
appropriate  to  these  discharges. 


UOLU 


Scattering  and  Transport  of  Resonance  Radiation  (Dr.  A.  V,  Phelps  and  Dr. 


J.  V.  Jenkins) 

During  the  period  of  this  report  we  have  completed  a  set  of  measure¬ 
ments  of  the  line  profiles  of  the  first  resonance  doublet  of  potassium  at 
moderately  high  densities.  We  have  been  especially  concerned  with  the 
accuracy  to  which  the  density  of  potassium  vapor  can  be  calculated  from 
the  equilibrium  vapor  pressure  data.  In  order  to  test  the  equilibrium 
data,  measurements  have  been  made  of  the  Intensity  of  radiation  reflected 
from  the  IntsrJace  between  the  potassium  vapor  and  the  glass  window.  The 
magnitude  of  the  reflection  coefficient  In  the  wings  of  the  profile  Is  a 
measure  of  the  product  of  the  ground  state  atom  density  and  the  oscillator 
strength  for  the  resonance  line.  Using  other  measurements  of  the  oscilla¬ 
tor  strength  one  can  determine  the  ground  state  density.  Alternately,  one 
can  use  the  equlllbr.'um  vapor  pressure  data  and  oscillator  strength  data 
to  predict  the  reflection  coefficient.  Figure  3  shows  such  a  comparison 
for  the  7644  A  line  of  potassium.  Because  of  some  discoloration  of  the 
windows  as  the  result  of  reactions  between  the  potassium  and  the  glass 
(Corning  1720) ,  there  Is  some  uncertainty  In  the  dielectric  coefficient 
of  the  glass.  We  are  therefore  repeating  the  measurements  with  an  optical 
arrangement  which  allows  much  more  rapid  collection  and  Intercomparison  of 


data. 
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.  Experimental  and  theoretical  specular  reflection  profiles  for 
6944  A  line  of  potassium. 
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Radlative  Hydrodynamics  (Dr.  D.  G.  Hummer,  Dr.  J.  Cassinelli,  Dr.  J. 

Castor  and  Mr.  P.  Kunasz) 

Line  formation  in  spherical  p;eomfc.:ry.  The  important  problem  of  spec¬ 
tral  line  formation  in  a  spherical  emitting  cloud  is  complicated  by  the 
existence  of  three  characteristic  scale  lengths,  which  may  be  taken  as 
line  center  total  optical  thickness,  the  thermalization  length  and  the 
radius  of  curvature.  In  more  conventional  plane  parallel  configuration 
only  the  first  two  are  encountered,  and  the  situations  occurring  when  one 
or  the  other  is  the  larger  are  well  understood.  However,  for  the  spherical 
case  very  little  is  understood  about  the  regime  defined  by  a  specific  ordering 
of  these  characteristic  lengths. 

As  an  aid  in  understanding  the  line  formation  process  in  spherical 
configurations  as  well  as  to  provide  the  capability  for  accurately  comput¬ 
ing  the  line  radiation  field  in  these  geometries,  we  have  generalized  the 

non-linear  iteration  procedure  reported  earlier  (the  so-called  variable 

1  2 

Eddington  factor  method)  for  monochromatic  transfer  ’  to  allow  for  non¬ 
conservative  transfer  in  a  gas  in  which  the  opacity  varies  with  frequency. 

The  coupling  of  radiation  of  different  frequencies  because  of  non-coherent 
scattering  is  taken  intcj  account.  Coupled  differential  equations  for  the 
monochromatic  mean  intensity,  containing  the  Eddington  factors  have  been 
derived  and  a  suitable  difference  scheme  has  been  devised.  A  code  embodying 
this  analysis  has  been  written  and  is  now  being  tested.  Some  trouble  has 
been  encountered  in  situations  in  which  the  optical  depth  is  very  small  and 
the  ratio  of  outer  to  inner  radii  is  very  large.  The  solution  of  the  dif¬ 
ficulty  is  now  being  pursued. 

In  conjunction  with  the  studies  on  line  formation,  the  subroutine 
FORMAL,  which  computes  the  formal  solution  of  any  spherical  transfer  pto- 
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blem  (given  the  opacity  and  source  function  as  functions  of  radius,  this 
subroutine  evaluates  the  radiation  field  and  its  first  three  moments)  has 
been  completely  rewritten  for  increased  accuracy  and  speed,  and  can  now 
operate  with  much  less  storage,  because  of  improved  logic  and  a  more  ac¬ 
curate  difference  scheme.  An  account  of  this  subroutine  is  being  prepared 
for  publication  in  Computer  Physics  Communications. 

Radiative  transfer  in  a  doublet  in  the  presence  of  excited  atom  dif¬ 
fusion.  As  a  further  extension  of  our  study  of  radiative  transfer  phen¬ 
omena  in  laboratory  situations,  we  have  generalized  our  previous  analysis 
of  radiative  transfer  in  a  spectral  line  in  situations  where  the  diffusion 
of  the  excited  atoms  is  Important  to  the  case  of  a  spectral  doublet.  We 
had  earlier  also  solved  the  doublet  transfer  problem  without  excited-state 
diffusion.  Because  of  the  very  great  complexity  of  the  present  problem, 
the  analysis  will  not  be  coded  until  a  clear  indication  of  the  need  for 
these  results  is  forthcoming. 

The  results  obtained  from  earlier  work  on  these  problems  is  now  being 
prepared  for  publication  and  the  possibility  of  publishing  a  set  of  com¬ 
puter  codes  for  laboratory  radiative  transfer  is  being  investigated. 

Radiation  hydrodynamics  jn  spherical  geometries.  The  transfer  of 
radiation  in  a  spherically  symmetric  moving  medium  has  been  considered  in¬ 
cluding  terms  of  order  v/c,  where  v  is  of  the  order  of  the  flow  velocity, 
using  a  formulation  in  which  the  radiation  quantities  are  those  in  the 
fluid  frame.  The  lowest  frequency-integrated  movement  equation  is  equiva¬ 
lent  to  the  first  law  of  thermodynajnlcs  for  the  radiation  and  its  use  in¬ 
sures  that  the  internal  energy  and  radiation  pressure  is  correctly  accounted 
for  in  the  conservation  of  energy  for  the  whole  system.  The  fluid  equations 
including  the  effects  of  matter-radiation  coupling  and  total  energy  con- 


servation  have  been  derived  in  a  form  which  shows  that  the  difference  be¬ 
tween  the  fluid  frame  and  stationary  frame  values  of  the  radiative  flux 
can  be  regarded  as  the  convective  flux  of  radiation  enthalpy.  This  work 
corrects  some  errors  in  earlier  treatments  of  this  problem.  Some  novel 
approximations  that  could  be  applied  to  various  fjims  of  the  transfer  equa 
tion  to  make  them  computationally  tractable  have  been  discussed  and  are 
currently  being  investigated,  including  further  generalizations  of  the 
variable  Eddington-f actor  method. 

References 

1.  D.  G.  Hummer  and  G.  B.  Rybicki,  Monthly  Notices  Roy.  Astron.  Soc.  152, 
1  (1971). 

2.  J.  P.  Cassinelli  and  D.  G.  Hummer,  Monthly  Notices  Roy.  Astron.  Soc., 
153,  9  (1971). 
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Plasma  Statistics  (Dr.  W.  E.  Brlttin  .and  Dr.  W.  R.  Chappell) 

Dr.  Chappell  is  presently  developing  a  more  general  kinetic  theory 
of  Stark  broadening  of  spectral  lines  in  order  to  refine  and  extend  the 
results  ob taint'd  previously  in  collaboration  with  Dr.  Cooper  of  JILA  and 
Dr.  Smith  of  NBS.'*’  This  work  involves  developing  a  kinetic  equation  for 
the  dipole  density  autocorrelation  function  in  a  manner  similar  to  that 
used  for  neutral  gas  broadening .  This  work  should  give  a  better  under¬ 
standing  of  the  effects  of  charged  particle  correlations.  A  paper  with 
R.  H.  Williams  was  recently  published  in  Physics  of  Fluids.^ 

References 

1.  W.  R.  Chappell,  J.  Cooper,  and  E.  W.  Smith,  "Non-thermal  Effects  in 
Stark  Broadening,"  J.  Quant.  Spectrosc.  Radiative  Transfer  1^,  1195 
(1970) . 

2.  W.  R.  Chappell,  J.  Cooper,  E.  W.  Smith,  and  T.  Dillon,  "A  Kinetic 
Theory  of  Spectral  Line  Shapes,"  J.  Stat.  Phys.  401  (1971). 

3.  W.  R.  Chappell  and  R.  H.  Williams,  "Microscopic  Theory  of  Conducti¬ 
vity  of  a  Weakly  Ionized  Plasma,"  Phys.  Fluids  1^,  1059  (1972). 
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‘II.  -TRANSPORT  COEFFICIENTS  AND  REACT lON  ilATES  ] 

I 

I 

’  Program  Leader  --Dr.  S.  J.  Smi,th  i 


'  '  JP' 


The  prediction  of  the  performance  characteristics  of  ionize, d  gas 

t 

devices,  such  as  the  efficiency,  power 'output  and  stability  of  high  power 
lasers,  requires  that  one  be  able  to  calculate  Accurate  values  of  elec¬ 
tron  transport  and  reaction  rate  coef f icieats  for  electrons  and  ions . '  In 
order  to  meet  some  of  these  requirements  we  afe  engaged  in;  q)  a  determin- 

t  ■  ' 

ation  of  sets  of  electron  collision  cto^s  sections  consistent  with  measured 

\ 

electron  transport  and  ionization  coefficients  and  with  recent  electron 

)  ' 

1  j 

beam  measurements  of  excitation  cross  sections  fot  the  upper  excited  states 

i  ,  , 

of  gases  found  in  molecular  lasers;  b)  measurement  of  cross  sections  for 

\ 

molecular  dissociation;  c)  evaluation  and  compilation  of  low  energy  electron, 

j 

photon,  and  ion  collision  cross  section  'date;  and  d)  determination  of  nega- 

!  '  '  !  ^  ' 
tlve  ion  characteristics,  such  as  single  ahd  multiple  pliotodetachment 'cross  ; 

i  j 

sections,  negative’ ion  stability  and  low  electron  energy  scattering  para¬ 


meters  . 


I  • 


t 


i 


! 


1 
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Electron  Transport  and  Ionization  Coefficients  (Dr.  G.  E.  Chamberlain, 

Dr.  L.  J  Kieffer  and  Dr.  A.  V.  Phelps) 

lor  the  past  six  mrnihf  our  efforts  have  been  concentrated  on  the 
development  of  criteria  to  be  used  in  the  selection  of  a  critical  set  of 
experimental  electron  transport  coefficients.  Emphasis  is  being  placed 
on  data  for  N2  and  on  high  values  of  E/N  for  which  electronic  excitation 
and  ionizaticn  are  important.  (E  =  electric  field,  N  *  gas  number  den¬ 
sity.)  The  Information  Center  bibliography  and  compilation  of  swarm  data 
being  prepared  by  J.  Dutton  has  helped  extensively.  Criteria  for  selec¬ 
tion  of  drift  velocity  data  W  are  currently  being  summarized.  The  drift 
velocity  work  will  be  followed  by  the  analysis  and  summary  of  criteria 
for  selection  of  data  on  D/y  (ratio  of  diffusion  to  mobility  coefficients), 
and  (ionization  coefficient) . 

Our  review  of  the  literature  shows  that  one  now  has  a  rather  high 
degree  of  understanding  of  errors  in  electron  drift  velocity  measurements 
caused  by  the  diffusion  of  electrons  along  density  gradients  resulting 
from  absorbing  boundaries,  ionization  growth  and  electron  injection  and 
observation  techniques.  Unfortunately,  only  a  few  experimeiitalists  have 
applied  these  corrections  cvr  have  given  sufficient  information  regarding 
experimental  conditions  so  that  we  can  apply  the  corrections.  It  appears 
that  the  situation  is  somewhat  less  satisfactory  with  regard  to  the  ef¬ 
fects  of  space  charge,  although  our  investigation  is  continuing.  It 
should  be  noted  that  the  careful  formulation  of  all  of  the  sources  of 
error  in  transport  coefficient.,  determinations  is  not  available  at  present 
so  that  this  effort  will  aid  others  in  making  more  reliable  measurements 
as  well  as  aid  us  in  selecting  data  for  use  in  our  cross  section  analysis. 

A  second  activity  during  this  reporting  period  has  been  inter :.ra- 
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parison  of  excitation  and  transport  coefficients  calculated  at  various 
laboratories  for  electrons  in  gas  mixtures  typical  of  electron  beam  sus" 
tained  lasers,  i.e.,  He:N2:C02  =  3;2:1.  The  objectives  of  such  a  com¬ 
parison  are  to  determine  whether  the  various  methods  for  solving  the  re¬ 
levant  Boltzmann  equation  are  equally  accurate,  to  insure  that  the  var¬ 
ious  groups  are  using  the  best  available  sets  of  electron  scattering  and 
excitation  cross  sections  and  to  establish  communication  with  the  var¬ 
ious  groups  so  that  improvements  in  cross  sections  or  computational  tech¬ 
niques  can  be  distributed  rapidly.  Inter comparisons  of  the  results  of 
our  calculations  with  those  of  three  other  laboratories  showed  one  set 
of  results  to  be  in  very  good  agreement  with  ours,  one  with  up  to  20%  too 
large  a  rate  of  vibrational  excitation  and  one  with  as  much  as  10%  toe 
low  a  rate  of  vibrational  excitation.  Sources  of  the  low  rate  of  excita¬ 
tion  are  being  investigated  through  further  calculations  by  the  laboratory 


involved. 
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JILA  Infoimatlon  Analysis  Center  (Dr.  L.  J.  Kieffer) 

The  Information  Center  prepared  a  report  on  the  ARPA-supported 
Workshop  "Dissociative  Excitation  of  Simple  Molecules"  which  was  held 
16-17  March,  1972.  The  report.  Issued  as  JILA  Information  Center  Report 
No.  12,  is  available  from  the  Defense  Documentation  Center,  Cameron  Sta¬ 
tion,  Alexandria,  Virginia  22314.  The  order  number  for  this  report  is 
AD745417. 

In  addition  to  abstracts  of  most  of  the  invited  papers,  the  report 
contains  an  extensive  compilation  of  cross  section  data  for  dissociative 
excitation  of  atmospheric  molecules  by  electron  impact  and  a  discussion 
of  the  molecular  branching  ratio  method  for  calibration  of  optical  sys¬ 
tems  in  the  vacuum  ultraviolet. 

Work  continues  on  the  publications  mentioned  in  the  last  semi-annual 
report.  The  publications  will  be  reported  when  available. 

Assistance  continues  to  be  given  to  Dr.  A.  V.  Phelps  in  his  attempt 
to  obtain  a  "best  set"  of  N2  cross  sections  for  use  in  laser  modeling  and 
modeling  of  other  devices  involving  low  pressure  ionized  gases . 
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Moieculag  Dissociation  Processes  (Dr.  L.  J.  Kieffer) 

During  the  month  of  July,  Professor  R.  J.  Van  Biiant,  University  of 
Virginia  (Visiting  Scientist  at  JILA) ,  took  extensive  data  on  two  dis¬ 
sociative  processes : 

(1)  e  +  O2  +  0  +  2e 

(2)  e  +  ©2  -*■  0  +  O"*"  +  e 

The  angular  distributions  of  the  ion  fragments  as  a  function  of  both 
the  ion  energy  and  the  electron  energy  for  both  these  processes  were  ob¬ 
served.  Also  observed  was  the  energy  distribution  of  the  0  ions  for  pro¬ 
cess  (2)  .(pair  formation)  as  a  function  of  the  electron  impact  energy. 

The  results  of  these  experiments  and  data  taken  previously  will  be  pre¬ 
pared  for  publication  within  the  next  six  months. 


Negative  Ion  Kinetics  (Dr.  W.  C.  Lineberger,  Dr.  J.  Hall,  Dr.  H.  Hotop, 

Mr.  R.  A.  Bennett,  Mr.  T.  A.  Patterson  and  Mr.  R.  A.  Beyer) 

We  have  continued  the  studies  of  metallic  negative  ions  produced  in 
our  sputtering  ion  source.  Photo  detachment  of  Au“  and  Pt“  has  been 
studied  using  a  tunable  dye  laser  with  a  1-2  A  bandwidth.  The  threshold 
behavior  of  the  photodetachu  »nt  cross  section  is  that  for  outgo:^Tig  p-wave 

electrons,  and  agrees  with  the  Wlgner  threshold  law,  a  «  k  over  at  least 

50  meV  above  threshold.  These  data  show  that  there  are  at  least  2  states 
-  9  2  2 

of  Pt  ,  and  thus  that  the  d  s  (  D)  configuration  is  stable.  Hartree- 

10  2 

Fock  estimates  by  Weber  at  Northwestern  Indicate  that  the  d  s(  S)  con¬ 
figuration  of  Pt  is  probably  unstable.  The  results  are  that  EA(Pt) 

2.128  ±  .002  eV  and  that  EA(Au)  »  2.3086  +  .0007  eV. 

In  a  separate  experiment  we  have  attached  our  sputtering  source  to 
the  laser-photoelectron  spectroscopy  apparatus  of  Dr.  J.  Hall,  in  which  one 
measures  the  energy  spectrum  of  electrons  detached  from  a  negative  ion  beam 
by  an  argon  ion  laser.  The  electron  affinities  of  Ag~  and  Cu  were  measured 
relative  to  0  and  OH  respectively;  the  Cu~  detached  electron  energy 
spectrum  is  shown  in  Fig.  4.  Our  results  are  that  EA(Ag)  =  1.303  ±  .008 
eV  and  EA(Cu)  =  1.226  ±  .008  eV.  All  of  these  electron  affinitie,;:;  are 
~  0.5  eV  less  than  surface  ionization  measurements,  and  the  EA(Cu)  is 
~  0.6  eV  less  than  estimates  based  on  accurate  Hartree-Fock  calculations 
with  an  empirical  electron  correlation  correction  added.  These  data,  which 
represent  the  first  accurate  determinations  of  electron  affinities  of  tran¬ 
sition  and  noble  metals,  make  it  clear  that  many  present  estimates  of  elec¬ 
tron  affinities  are  Incorrect.  We  hope  to  measure  a  number  of  additional 
electron  affinities  of  metals  in  the  near  future. 
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We  have  also  used  the  tunable  dye  laser  to  investigate  photodetach¬ 
ment  of  OH  and  OD  in  the  range  7000A  -  6400A.  A  detailed  analysis  of 

the  observed  cross  sections  is  being  made  on  the  basis  of  the  known  term 

2 

values  for  OH  and  OD  (X  and  the  experimental  result  of  Hall,  Celotta 

and  Bennett  that  »0D  )  =  Yg(OH,OD)  to  within  i.OOlA.  The  electron 

affinities  are  found  to  be  EA(OH)  *  1.826  ±  .002  eV,  and  EA  CD)  =  1.823  ± 
.002  eV.  These  values  are  in  agreement  with  Branscomb's  result  of  1.83  ± 

.04  eV  as  extracted  from  a  low  resolution  experiment. 


ELECTRON  BINDING  ENERGY,  eV 


fig.  4.  Electron^energy  spectrum  for  electrons  detached  from  Cu  ions 
by  4880  A  photons. 
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III.’~>ENERGY  Lt)SS  AND  RADIATION^PRODUCTION  BY  ELECTRON  COLLISIONS  ' 

i 

Program  Leader  -  Dr.  E.  C.  Beaty 


Interaction  of  Electrons  with  Ions  (Dr.  G.  H.  Dunn,  Dr.  F.  L.  Walls,  Dr. 

D.  Crandall  and  Mr.  P.  0.  Taylor) 

The  interactions  of  electrons  with  ions  can  be  very  Important  in  par- 
tially  i07il2ed  gases  because  of  the  very  large  cross  sections  for  scatter- 
Ing,  exv-itatlon  and  dissociative  recombination.  We  are  engaged  in  measure¬ 
ments  of  excitation  cross  sections  and  recombination  coefficients  for  ions 
of  atmospheric  interest  and  of  importance  for  testing  theory.  The  energy 
dependence  of  electron-molecular  ion  recombination  coefficients  has  taken 
on  added  importance  with  the  development  of  electron  beam  excited  and  sus¬ 
tained  lasers. 

The  cross  sections  for  electron-impact  excitation  of  N2"*’,  Ca"*',  and 
Ba  are  measured  using  crossed  beams  of  variable  energy  electrons  and  the 
respective  ions.  The  resultant  resonance  radiation  is  observed  along  the 
third  perpendicular  direction.  Preliminary  results  were  reported  earlier. 
During  this  report  period  absolute  radiometry  was  completed  for  the  mea¬ 
surements  cited  above.  Final  cross  section  measurements  were  made  on 
Ca  from  2  eV  to  700  eV,  and  corrections  to  the  data  were  carefully  ana¬ 
lyzed  and  applied.  Figure  5  shows  a  plot  of  cross  section  versus  elec¬ 
tron  energy  for  both  the  H  (4^P^^2  ^  ^  ^^^1/2^ 

emissions  at  3968  and  3934  A  respectively.  The  measurements  are  abso¬ 
lute,  and  the  flags  represent  98%  confidence  limits  (about  3  standard  de¬ 
viations)  .  Points  plotted  with  a  B  represent  the  three-state  close  coup¬ 
ling  calculation  of  Burke  and  Moores.^  Points  plotted  with  S  represent 
a  Coulomb  distorted  wave  calculation  of  Burgess  and  Sheory.  The  measure- 
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mants  indicate  that  further  theory  and/or  cociputation  is  needed  before 
such  calculations  can  be  trusted  to  within  30%  in  the  threshold  region 
where  such  cross  sections  are  of  value  in  understanding  various  plasmas. 

Figure  6  shows  a  plot  of  the  measured  polarization  for  the  K  line 
as  a  function  of  electron  energy.  The  points  plotted  with  an  S  are  due 
to  a  calculation  by  Saraph  using  Burke  and  Moores  close-coupling  matrix 
elements.  Here  the  flags  represent  one  standard  deviation  of  statistical 
fluctuations  combined  with  systematic  uncertainties.  Again  the  measure¬ 
ments  indicate  a  need  for  further  theoretical  development.  Work  is  con¬ 
tinuing  on  the  Ba^  and  N2^  measurements . 

The  experiment  for  measuring  the  energy  dependence  of  electron 
molecular  ion  recombination  coefficients  employs  trapping  of  ions  in  a 
dc  electric  and  magnetic  field  configuration  to  investigate  collisions  of 
these  ions  with  other  particles.  In  particular  dissociative  recombination 
of  these  ions  with  monoenergetic  electrons  will  be  investigated  over  an 
electron  energy  range  from  0  to  several  eV.  This  report  period  was  spent 
in  apparatus  improvements  looking  forward  to  serious  data  taking.  A  new 
trap  was  fabricated  from  molybdenum,  and  its  closer  shape  tolerance  and 
low,  uniform  surface  potentials  should  result  in  a  device  with  narrower 
line  widths.  An  interface  was  designed  for  running  the  experiment  with 
a  small  computer,  and  programs  were  written  to  do  so.  Extra  pumping  was 
added  to  the  system  both  in  the  form  of  liquid  air  trapping  and  Ti  sub¬ 
limation  pumping.  This  resulted  in  a  base  pressure  of  2  x  10  torr,  in¬ 
creased  ion  residence  times ,  and  decreased  ion-molecule  reactions  to  an 
apparently  acceptable  level.  Apparatus  was  assembled  for  testing  the 
Channeltron  bundle  detector  designed  to  measure  the  spatial  distribution 
of  ions  in  the  trap. 


ELECTRON  ENERGY  (eV) 

Fig.  5.  Absolute  cross  sections  for  exciting  the  k  and  H  Wssions. 


Present 

Experiment 


K  line  a(V^  -  -  av) 

K  line  a(V  -  <}>) 

Vf 

H  line  a(Vg  -  -  AV) 


IB  Close  coupling  calculation  6f  Burke  and  Moores 
(Ref.  1). 

S  Coulomb  distorted  wave  calculation  of  Burgess  and 
Sheory  (Ref.  2).  ' 

pressed  at  the  98/.  confidence  level.  Above  60  eV  for  the  K  line  an  addi 
f  uncertainty  has  been  allowed  for  the  possible  presence  o 

error  ba?  ^  ^he  inner 

lute  opti^/rairratio:^ 


Polarization  of  the  K  emission.  The  inset  shows  the  low  energy  region  on  an  expanded  scale. 
Flags  represent  one  S.D.  of  statistical  fluctuations  combined  with  systematic  errors. 

X  -  Experimental  data  taken  in  1972;  0  -  Experimental  data  taken  in  1970;  -S-  Theoretical 

polarization  from  Saraph  (Ref.  3)  based  on  the  close  coupling  calculation  of  Burke  and  Moores 


Theoretical  Atomic  Physics  (Dr.  D.  G.  Hummer,  Dr.  J.  I.  Castor  and  Dr. 

D.  W.  Norcross) . 

The  main  emphasis  in  this  period  has  been  on  t  calculation  of  elec¬ 
tron-excitation  cross  sections  for  CIII  and  NIV  considering  all  transitions 
among  the  lowest  14  terms.  Eleven  configurations  are  included.  Soiie  pro¬ 
gress  has  been  made  on  the  lithium-like  ions  CIV,  KV  and  OVI  using  9-term 
models.  These  cross  sections  are  being  calculated  in  the  unitarized  dis¬ 
torted-wave  calculation.  Much  of  the  effort  in  this  period  has  gone  i.nto 
improving  the  efficiency  of  the  programs;  as  well,  a  considerable  number 
of  R-matrix  elements  have  been  calculated,  although  no  cross  sections  have 
been  finished  to  date. 

We  now  have  the  capability  to  calculate  resonance  contributions  to 
these  cross  sections  which  are  ignored  in  the  usual  distorted  wave  proce¬ 
dure.  As  sliown  by  Seaton,^  multichannel  quantum-defect  theory  allows  one 
to  use  certain  smoothly  varying  energy-dependent  matrices,  which  can  be 
computed  in  the  distorted-wave  approximation,  to  compute  the  resonance 
structure  in  electron-scattering  cross  sections.  The  code  inccrporating 
this  analysis  was  written  in  London  by  a  member  of  Professor  Seaton's 
group  and  as  a  result  of  Hummer's  recent  visit  to  London  and  Seaton's 
visit  to  JILA,  this  facility  is  now  operational  here.  It  reproduces  the 
results  of  coupled-equations  calculations  using  about  1/100  of  the  compu¬ 
ter  time  used  by  previously  available  techniques.  The  handling  and  final 
processing  of  the  vast  quantities  of  intermediate  data  involved  now  poses 
a  serious  problem  and  considerable  effort  is  being  spent  on  the  design  of 
an  automatic  data-handling  procedure  that  is  convenient  yet  reliable  des¬ 
pite  computer  unreliability. 

Reference 

1.  M.  J.  Seaton,  J.  Phys.  B  (Proc.  Roy.  Soc.)  5  (1969). 
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Atomlc  Interaction  Theory  (Dr.  S.  Geltman  and  Dr.  D.  W,  Norcross) 

A  general  computer  program  has  been  prepared  to  calculate  free-free 
absorption  coefficients  (inverse  bremsstrahlung)  and  emissivities  for 
neutral  atoms.  The  continuum  atomic  wave  functions  are  based  on  the 
Robinson-Geltman  potential^ 


n-u  • 

(t/  +  r2)^ 


U) 


where  U„_  is  the  Hartree-Fockr Slater  potential  as  evaluated  by  the  me- 
Ho 

2 

thods  of  Herman  and  Skillman,  the  second  tern  removes  the  Coulomb  tail 

as  «ieen  by  an  electron  within  the  neutral  atom,  apd  the  third  term  is  the 

polarization  potential.  The  best  available  polarizabilities  are  taken 

3 

from  a  recent  complete  tabulation,  and  the  r^  and  r^  parameters  are 
adjusted  to  give  the  best  fit  with  elastic  scattering  data  or  electron 
affinities.  All  dipole  matrix  elements  are  evaluated  numerically  and 
enough  partial  waves  are  kept  to  insure  convergence  of  results.  Results 
are  obtained  for  the  atmospheric  atoms  N  and  0,  as  well  as  C  and  all 
the  rare  gases.  The  ranges  of  temperature  and  wavelength  covered  are 
T  ®  500°  to  20000°K,  X  =  ,5  to  20  microns.  Comparison  with  the  shock 

4 

tube  studies  of  Taylor  and  Caledonia  shows  good  agreement  for  the  case 
of  Ne,  while  the  shock  tu'be  results  are  about  a  factor  of  three  higher 
than  the  present  calculation  for  N,  0,  A,  and  Xe.  We  estimate  the  cal¬ 
culation  to  be  good  to  about  30%  on  the  basis  of  theoretical  consistency 
and  our  fits  to  the  scattering  data.  We  plan  to  evaluate  the  free-free 
absorption  coefficients  of  the  metallic  species  Mg,  Ae,  Si,  Fe,  and  Ce. 
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Table  IJ  Free-Free  Atomic  Absorption  i 
Coefficients  in  Units  of  cm^. 


Atom  and  Temperature. 

X  in  Microns 

Present  Calculation 

Taylor  &  Caledonia 

2.0 

,  ’  8.86  (-39) 

3.20  (-38) 

Nitrogen  9175  K 

3.5 

3.97 ’(-38)  ,  .  ' 

1.30  (-37) 

5.0 

1 

1.08  (-37) 

'3.00  (-37) 

o 

r  2.0  . 

6^82  (-39) 

1  1 

2.90  (-38)  ' 

Oxygen  9700  K 

3.5  : 

-  L 

3.12  (-38) 

,1.00  (-37) 

5.0 

8.53  (-38)  , 

1 

2.60  (-37) 

2.0 

5.40  (-39)  1 

5.80  (-39) 

Neon  12600  K 

t 

3.5 

2.59  (-38) 

2.70  (-38) 

5.0 

7.21  (-38i) 

7.2q  (-38) 

o 

2.0 

8.34  (-39) 

3.80  (-38) 

Argon  9900  K 

3.5 

3.81  (-38 j 

'  1.50  (-37)  ' 

5.0 

1^04  (-37)  , 

1 

3.60  (-37) 

2.0 

1.51  (-38) 

3.00  (-38) 

Xenon  8200  K 

3.5 

6.50  (-38)  '  ■ 

9.60  (-38) 

5.0 

I 

1.74  (-37) 

2.00  (-37) 

i  [ 
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Electron  scattering  and  photolonlzatlon.  A  psper  on  the  close— coup¬ 
ling  calculations  of  cross  sections  for  low-energy  (0-5  eV)  scattering 
of  electrons  by  sodium  has  been  accepted  for  publication.^  Additional 
experimental  verification  of  these  results  has  recently  been  obtained. 

,  These  results  are  confined  to  elastic  scattering  from,  and  excitation  of, 
the  ground  (3s)  state.  We  are  now  proceeding  to  extract  cross  sections 
for  elastic  scatter'‘.ng  from,  excitation  of,  and  transitions  between  the 
substates  3p(j=l/2)  and  3p(j=3/2)  of,  the  first  excited  state.  Cross 
sections  of  the  latter  type  are  of  particular  Importance  in  the  develop¬ 
ing  field  of  infrared  astronomy. 

1  Reliable  estimates  of  the  photodetachment  crb3s  sections  and  elec¬ 
tron  affinities  of  the  negative  ions  of  the  alkali  metals  are  required 
for  the  interpretation  of  the  properties  of  low-temperature'  plasmas.  No 
reliable  experimental  data  are  available  on  the  electron  affinities  of 
the  alkali  metals  and  none  whatsoever  on  the  photodetachment  cross  sec¬ 
tions.  Encouraged  by  the  above  results  for  sodium,  photodetachment  cross 
sections  of  the  lithium  and  sodium  negative  ions  have  been  computed  for 
ejected  electron  energies  from  threshold  to  5  eV  (X  ~  20,000  to  2000  A), 
using  configuration-interaction  (Cl)  wave  functions^  for  the  initial 
states  and  solutions  of  the  close-couplir»g  equations  for  the  final  states. 
The  results^  satisfy  two  "necessary"  conditions  quite  well,  and  are  signi¬ 
ficantly  different  from  those  of  much  less  sophisticated  calculations.  A 
sharp  resonant  feature  occurs  in  both  cross  sections  at  wavelengths  cor¬ 
responding  to  transitions  to  the  first  excited  state  of  the  neutral  atom. 
This  feature  is  so  pronounced  that  it  could  be  used  for  a  precise  deter- 
mination  of  the  electron  affinities  of  lithium  and  sodium  using  present 
dye  laser  techniques.  Equally  reliable  Cl  wave  functions  are  available 
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for  the  negative  ion  of  potassium,^  and  these  calculation?  will  be  extended 
to  photodetachnent  of  this  species  as  soon  as  a  suitable  target  model  for 
the  electron-potassium  scattering  calculation  is  obtained. 

7 

A  study  of  photoabsorption  by  cesium  has  been  completed.  The 

results  confirm  the  importance  of  the  spin-orbit  and  core-polarization 

effects,  and  demonstrate  that  relativistic  and  aon-adiabatic  corrections 

0 

omitted  from  previous  work  are  not  negligible.  They  also  lend  no  support 

9 

to  a  recent  suggestion  that  there  is  a  pole  in  th^'.  doublet  line  strength 
ratio  for  principal  series  transitions  at  principal  quantum  number  n  ~  10 
15,  and  suggest  that  the  accepted  photoionization  cross  section^®  is  too 
large  by  about  a  factor  of  two.  This  work  will  be  the  basis  of  a  close¬ 
coupling  calculation  of  low-energy  electron  scattering  by  cesium,  in  an 

effort  to  resolve  outstanding  discrepancies  between  experimental  and  theo- 
11  12 

retical  result®.  *  A  similar  study  of  photoabsorption  by  potassium 
will  be  undertaken  for  reasons  discussed  above. 
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lonlzlng  Collisions  (Dr.  E.  C.  Beaty  and  Dr.  J.  Hauer) 

Work  was  restarted  on  the  project  to  study  the  energy  and  angular 
distributions  of  secondary  electrons  with  the  arrival  of  Dr.  John  Mauer. 
Some  pieces  of  equipment  failed  during  the  shutdown,  which  combined  with 
the  new  personnel  has  required  redoing  many  of  the  equipment  verification 
checks.  This  restarting  procedure  is  nearly  finished.  Some  extensions 
of  the  earlier  data  are  planned.  A  new  system  for  mounting  the  electron 
gun  and  a  second  detector  have  been  designed  and  built.  These  new  parts 
cannot  be  Installed  in  the  vacuum  system  until  the  present  sequence  of 
measurements  is  completed. 
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IV.  ^GENERATION  AND  ABSORPTION  OF  .^lATION  , 


Program  Leader  -  D»-.  A.  C.  Gallagher 


Alkali-Inert  Gas  Line  Broadening  and  Intermolecular  Potentials  (Dr.  A. 

C.  Gallagher  and  Mr.  D.  L.  Drummond) 

The  measurements  of  rubidium- inert  gas  spectra  have  been  completed,, 
and  are  being  reduced  to  obtain  the  interaction  potentials.  In  addition, 
undulations  in  these  continuous  spectra  have  been  detected  as  shown  by 
Fig.  7.  These  are  attributed  to  transitions  from  bound  vibrational  states 
in  the  excited  electronic  state.  Computations  of  bound -bound  and  bound- 
free  Franck-Condon  factors  are  now  underway  to  utilize  these  data  as  a 
definite  test  of  the  adiabatic  potentials  Inferred  from  the  continuum 
data.  Since  these  maxima  in  the  continuum  spectra  Identify  vibrational 
states  of  the  excited  electronic  state  (the  pressure  dependence  of  the  am¬ 
plitudes  of  the  maxima)  can  be  used  to  test  models  for  three-body  mole¬ 
cular  formation. 
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Figure  7.  Emission  spectrum  of  the  rubidium  resonance  line  perturbed  by 
inert  gases.  The  undulations  are  attributed  to  bound  free 
transitions. 
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Plasma  Spectroscopy  (Dr.  J.  Cooper,  Mr.  J.  Bauer,  Mr.  S.,  A.  Freudenstein 
and  Mr.  J.  Slater) 

The  theoretical  work  on  the  redistribution  of  radiation  in  the  pre¬ 
sence  of  collisions  has  been  extended  to  take  into  account  the  effects  of 
a  magnetic  field. ^  The  results,  which  appear  complicated,  do  however, 
lend  themselves  to  simple  physical  interpretation  in  terns  of  magnetic 
field  depolarization  and  multiple  collision  rates. 

The  tables  on  hydrogen  line  broadening  have  been  accepted  for  publi- 
2 

cation,  and  recent  work  on  Stark  broadening  has  concentrated  on  inclu¬ 
sion  of  ion  dynamics.  Preliminary  results  give  excellent  agreement  with 
the  Hel  AA70  A  forbidden  line.  See  Fig.  8. 

3 

Proton  broadening  of  the  wings  of  Lyman  a  has  also  been  studied 
using  accurate  H2'*'  potentials.  Of  particular  interest  are  various  satel¬ 
lites  which  occur  on  the  red  wing  at  1233.5,  1240.5  and  1409&,  as  shown 
in  Fig.  9. 

Experimental  work  during  this  period  has  involved  rebuildii g  equip- 
ment  for  obtaining  high  densities  (c:  10  cm  ”)  of  H  for  absorption  ex¬ 
periments  . 
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Fig  9.  Absorp'*;ion  coefficient  a  for  red  wing  of  Lyman  a.  The  Holtzmark 
profile  is  shown. 
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Predlctlon  of  Laser  Gain  (Dr.  A.  V.  Phelps) 

The  theory  and  understanding  developed  during  the  study  of  alkali- 
inert  gas  broadening  and  intermolecular  potentials,  as  discussed  above, 
has  been  applied  to  the  prediction  of  the  gain  of  -'roposed  laser  systems 
utilizing  molecules  with  dissociating  ground  states,  e.g.,  Xe2  of  LiXe. 
Thus,  we  have  utilized  theoretically  or  experimentally  determined  inter¬ 
molecular  potential  energy  curves  for  the  ground  state  and  first  excited 
state  of  molecules  to  calculate  the  magnitude  and  spectral  variation  of 
absorption  and  stimulated  emission  coefficients.  These  calculations  show 

that  the  densities  of  ‘excited  molecules  required  for  lasers  of  reasonable 

15  3 

length  are  high,  e.g.,  of  the  order  of  10.  molecules/cm  ,  but  are  appar¬ 
ently  attainable  with  sufficiently  intense  excitation.  Figure  10  shows 
results  obtained  for  the  LlXe  system.  A  detailed  discussion  of  these  re¬ 
sults  has  been  prepared  and  will  be  sent  to  those  on  the  ARPA  distribution 
list  as  JILA  Report  110. 
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Fig.  10.  Stimulated  emission  and  fluorescence  in  xenon. 
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